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We show how to compute or at least to estimate vari- 
ous capacity-related quantities for Bosonic Gaussian channels. 
Among these are the coherent information, the entanglement 
assisted classical capacity, the one-shot classical capacity, and 
a new quantity involving the transpose operation, shown to be 
a general upper bound on the quantum capacity, even allowing 
for finite errors. All bounds are explicitly evaluated for the 
case of a one-mode channel with attenuation/amplification 
and classical noise. 



I. INTRODUCTION 

During the last years impressive progress was achieved 
in the understanding of the classical and quantum capac- 
ities of quantum communication channels (see, in partic- 
ular, the papers 0] - ||, where the reader can also find 
further references). It appears that a quantum channel 
is characterized by a whole variety of different capacities 
depending both on the kind of the information transmit- 
ted and the specific protocol used. 

Most of this literature studies the properties of sys- 
tems and channels described in finite dimensional Hilbert 
spaces. Recently, however, there has been a burst of in- 
terest (see e. g. 0) in a new kind of systems, sometimes 
called "continuous variable" quantum systems, whose ba- 
sic variables satisfy Heisenberg's Canonical Commuta- 
tion Relations (CCR). There are two reasons for this new 
interest. On the one hand, such systems play a central 
role in quantum optics, the canonical variables being the 
quadratures of the field. Therefore some of the current 
experimental realizations Q of quantum information pro- 
cessing are carried out in such systems. In particular, 
the Bosonic Gaussian channels studied in this paper can 
be seen as basic building blocks of quantum optical com- 
munication systems, allowing to build up complex opera- 
tions from "easy, linear" ones and a few basic "expensive, 
non-linear" operations, such as squeezers and parametric 
down converters. 

The other reason for the interest in these systems is 
that in spite of the infinite dimension of their underly- 
ing Hilbert spaces they can be handled with techniques 



from finite dimensional linear algebra, much in anal- 
ogy to the finite dimensional quantum systems on which 
the pioneering work on quantum information was done. 
Roughly speaking this analogy replaces the density ma- 
trix by the covariance matrix of a Gaussian state. Then 
operations like the diagonalization of density matrices, 
the Schmidt decomposition of pure states on compos- 
ite systems, the purification of mixed states, the com- 
putation of entropies, the partial transpose operation on 
states and channels, which are familiar from the usual 
finite dimensional setup, can be expressed once again by 
operations on finite dimensional matrices in the contin- 
uous variable case. The basic framework for doing all 
this is not new, and goes under heading "phase space 
quantum mechanics" or, in the quantum field theory and 
statistical mechanics communities, "quasi-free Bose sys- 
tems" H . Both authors of this paper have participated in 
the development of this subject a long time ago [|To|-p"2[ . 
In this paper, continuing |13| and H], we make further 
contributions to the study of information properties of 
linear Bosonic Gaussian channels. We focus on the as- 
pects essential for physical computations and leave aside 
a number of analytical subtleties related to infinite di- 
mensionality and unboundedness unavoidably arising in 
connection with Bosonic systems and Gaussian states. 

The paper is organized as follows. In the Section II we 
recapitulate some notions of capacity, which are currently 
under investigation in the literature, and what is known 
about them. Naturally this cannot be a full review, but 
will be limited to those quantities which we will evaluate 
or estimate in the subsequent sections. A new addition 
to the spectrum of capacity-like quantities is discussed 
in Subsection II. B: an upper bound on the quantum ca- 
pacity (even allowing finite errors), which is both simple 
to evaluate and remarkably close to maximized coherent 
information, a bound conjectured to be exact. In Sec- 
tion III we summarize the basic properties of Gaussian 
states. Although our main topic is channels, we need 
this to get an explicit handle on the purification opera- 
tion, which is needed to compute the entropy exchange, 
and hence all entropy based capacities. Bosonic Gaussian 
channels are studied in Section IV. Here we introduce the 
techniques for determining the capacity quantities intro- 
duced in Section I, deriving general formulas where pos- 
sible. In the final Section V we apply these techniques 
to the case of a single mode channel comprising attenu- 
ation/amplification and a classical noise. Some technical 
points are treated in the Appendices. 
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II. NOTIONS OF CAPACITY 

A. Basic entropy and information quantities 

Consider a general quantum system in a Hilbert space 
H = Hq. Its states are given by density operators p on 
Ji. A channel is a transformation p — > T[p] of quantum 
states of the system, which is given by a completely posi- 
tive, trace preserving map on trace class operators. This 
view of channels corresponds to the Schrodinger picture. 
The Heisenberg picture is given by the dual linear oper- 
ator X — > T*[X] on the observables X, which is defined 
by the relation 

TtT[p]X = Tr pT*[X], 

and has to be completely positive and unit preserving (cf. 

©)■ ' 

It can be shown (see e.g. flq]) that any channel in this 

sense arises from a unitary interaction U of the system 
with an environment described by another Hilbert space 
He which is initially in some state p E , 

T[p]=Tr E U(p®p E )U*, 

where Tt e denotes partial trace with respect to He, and 
vice versa. The representation is not unique, and the 
state pe can always be chosen pure, pe — |'0s)(V'-e|- The 
definition of the channel has obvious generalization to the 
case where input and output are described by different 
Hilbert spaces. 
Let us denote by 



H(p) = -Trp log p 



(2.1) 



the von Neumann entropy of a density operator p. We 
call p the input state, and T[p] the output state of the 
channel. There are three important entropy quantities 
related to the pair (p, T) , namely, the entropy of the in- 
put state H(p), the entropy of the output state H(T[p\), 
and the entropy exchange H(p,T). While the definition 
and the meaning of the first two entropies is clear, the 
third quantity is somewhat more sophisticated. To de- 
fine it, one introduces the reference system, described by 
the Hilbert space TLr : isomorphic to the Hilbert space 
TLq = TL of the initial system. Then according to |fl7| , 
H, there exists a purification of the state /?, i.e. a unit 
vector \tp) € Hq (g> TLr such that 

p = Tr R \ip){ip\. 

The entropy exchange is then defined as 

H(p,T) = H{(T®id)[\ip){i>\]), (2.2) 

that is, as the entropy of the output state of the dilated 
channel (Tig) id) applied to the input which is purification 
of the state p. Alternatively, 

H(p,T) = H(p' E ), 



where p' E = Te[p] is the final state of the environment, 
and the channel Te from TLq to TLe is defined as 

T E [p] = Ti Q U(p®p E )U*, 

provided the initial state pe of the environment is pure 
@, §• 

[From these three entropies one can construct several 
information quantities. In analogy with classical informa- 
tion theory, one can define quantum mutual information 
between the reference system R (which mirrors the input 
Q) and the output of the system Q' [[[7) , Q as 



I(p,T)=H(p' R )+H(p' Q )-H(p' RQ ) 
= H(p)+H(T[p])-H(p,T). 



(2.3) 



The quantity I(p, T) has a number of nice and "natural" 
properties, in particular, positivity, concavity with re- 
spect to the input state p and additivity for parallel chan- 
nels [||. Moreover, the maximum of I(p,T) with respect 
to p was argued recently to be equal to the entanglement- 
assisted classical capacity of the channel ||, |18|, namely, 
the classical capacity of the superdense coding protocol 
using the noisy channel T . It was shown that this max- 
imum is additive for parallel channels, the one-shot ex- 
pression thus giving the full (asymptotic) capacity. 

It would be natural to compare this quantity with 
the (unassisted) classical capacity C(T) (the definition 
of which is outlined in the next Subsection); however 
it is still not known whether this capacity is additive 
for parallel channels. This makes us focus on the one- 
shot expression, emerging from the coding theorem for 
classical-quantum channels Q 



Ci(T) 



(2.4) 



where the maximum is taken over all probability distri- 
butions {pi} and collections of density operators {pi} 
(possibly satisfying some additional input constraints). 
Ci(T) is equal to the capacity of T for classical infor- 
mation, if the coding is required to avoid entanglement 
between successive inputs to the channel. The full ca- 
pacity is then attained as the length n of the blocks, over 
which encoding may be entangled goes to infinity, i.e., 



C{T) = lim -d(T® n ). 

n — >oq ft 



(2.5) 



An important component of I(p, T) is the coherent in- 
formation 



J(p,T) = H(T[p})-H(p,T), 



(2.6) 



the maximum of which has been conjectured to be the 
(one-shot) quantum capacity of the channel T |ljj, ||. 
Its properties are not so nice. It can be negative, its 
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convexity properties with respect to p are not known, 
and its maximum was shown to be strictly superadditive 
for certain parallel channels |2(]], hence the conjectured 
full quantum capacity may be greater than the one-shot 
expression, in contrast to the case of the entanglement- 
assisted classical capacity. In this paper we shall also 
compare this expression with a new upper bound on the 
quantum capacity Q(T) (as introduced e.g. in the next 
Subsection). 

B. A general bound on quantum channel capacity 

In this Subsection we will establish a general estimate 
on the quantum channel capacity, which will then be eval- 
uated in the Gaussian case, and will be compared with 
the estimates of coherent information. Let us recall first 
a definition of the capacity Q(T) of a general channel T 
for quantum information. Intuitively, it is the number 
of qubits which can be faithfully transmitted per use of 
the channel with the best possible error correction. The 
standard of comparison is the ideal 1-qubit channel id2, 
where id„ denotes the identity map on the n x n-matrices. 
Then the quantum capacity Q(T) of a channel T (possi- 
bly between systems of different type) is defined as the 
supremum of all numbers c, which are "attainable rates" 
in the following sense: For any pair of sequences n a , m a 
with lim a (n a /m a ) = c we can find encoding operations 
E a and decoding operations D a such that 

\\idf n * -D a T® m «E a \\ ch ^0- 

Here || • jjcb is the so-called "norm of complete bound- 
edness" plfl , which is defined as the supremum with re- 
spect to n of the norms ||(T <g) id„)|] . It is equal to the 
"diamond metric" introduced in p2]. We use this norm 
because on the one hand, it leads to the same capacity as 
analogous definitions based on other error criteria (e.g., 
fidelities f|||) and, on the other hand, it has the best 
properties with respect to tensor products, which are our 
main concern. In particular, \\T ® S\\ c b = ||T|| c b ■ ||<S|| cb . 
Completely positive maps satisfy ||T|| c b = \\F\\, where F 
is the normalization operator determined by Tr(T[p]) = 
Tr(pF). In particular, ||T|| c b = 1 for any channel. We 
also note another kind of capacity, in which a much 
weaker requirement is made on the errors, namely 

||idf n ° -D a T® m *E a \\ ch <e<l (2.7) 

for all sufficiently large a, and some fixed e. We call the 
resulting capacity the e-quantum capacity, and denote 
it by Q e (T). Of course, Q(T) < Q £ {T), and by anal- 
ogy with the classical case (strong converse of Shannon's 
Coding Theorem) one would conjecture that equality al- 
ways holds. 

The unassisted classical capacity C(T) can be defined 
similarly with the sole difference that both the domain 
of encodings E and the range of decodings D should be 



restricted to the state space of the Abelian subalgebra 
of operators diagonalizable in a fixed orthonormal basis. 
In that case there is no need to use the cb-norm, as it 
coincides with the usual norm. According to recently 
pro ven strong converse to the quantum coding theorem 
@, H, C e (T) = C{T) where C e {T) is defined similarly 
to Q e {T). 

The criterion we will formulate makes essential use of 
the transpose operation, which we will denote by the 
same letter O in any system. For matrix algebras, 
can be taken as the usual transpose operation. How- 
ever, it makes no difference to our considerations, if any 
other anti-unitarily implemented symmetry (e.g. time- 
reversal) is chosen. In an abstract C*-algebra setting 6 
is best taken as the "op" operation, which maps every 
algebra to its "opposite" . This algebra has the same un- 
derlying vector space, but all products AB are replaced 
by their opposite BA. Obviously, a commutative alge- 
bra is the same as its opposite, so on classical systems 
O is the identity. Although the transpose maps density 
operators to density operators, it is not an admissible 
quantum channel, because positivity is lost, when cou- 
pling the operation with the identity transformation on 
other systems, i.e., O is not completely positive. A sim- 
ilar phenomenon happens for the norm of O: we have 
||0|| c b > 1 unless the system is classical. In fact, 

||e n || eb = n, (2.8) 

where ©„ denotes the transposition on the n x n-matrices 
pl[ . We note that since we do not distinguish the trans- 
pose on different systems in our notation, the observation 
that tensor products can be transposed factor by factor 
is expressed by the equation O = O (g) 9. Moreover, 
although for a channel T, the operator TO may fail to 
be completely positive, ©TO is again a channel, and, in 
particular, satisfies ||OTO|| c b = 1. 

The main result of this Subsection is the estimate 

Qe{T) < log||T6|| cb = Qe(T), (2.9) 

for any channel T. The proof is quite simple. Suppose 
n a /m a — > c < Q E (T), and encoding E a and decoding D a 
are as in the definition of Q £ (T). Then by Equation ( pT§| ) 
we have 

2 n a = || id |n a e|| cb < 

< ||(idf n ° -D a T® m °E a )Q\\ ch 

+ \\D a T® m "E a e\\ ch 

< ||e a »a||cb ||idf - D a T® m "E a \\ ch 

+\\D a (Te)® m «OE a e\\ cb 

<2"° £ + ||TO||™", 

where at the last inequality we have used that D a and 
QE a Q are channels, and that the cb-norm is exactly ten- 
sor multiplicative, so ||A® m || cb = ||^||" b - Hence, by tak- 
ing the logarithm and dividing by m Q , we get 
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log2 



log(l - e) 



<log||T6| 



cb' 



If we take base 2 logarithms, as is customary in informa- 
tion theory, we have log 2 = 1. Then in the last inequality 
we can go to th e limit a — > oo, obtaining c < Q@(T), and 
Equation (2.9) follows by taking the supremum over all 
attainable rates c. Note that base 2 logarithms are built 
into the above definition of capacity, because we are us- 
ing the ideal qubit channel as the standard of comparison. 
This amounts only to a change of units. If another base 
is chosen for logarithms is chosen, this should also be 
done con sistently in all entropy expressions, and Equa- 
tion (2.9) holds once again without additional constants. 

The upper bound Q&(T) computed in this way has 
some remarkable properties, which make it a capacity- 
like quantity in its own right. For example, it is exactly 
additive: 



Qe(S®T) = Qe{S) + Q e (T), 



(2.10) 



for any pair S, T of channels, and satisfies the "bottleneck 
inequality" Q e (ST) < mm{Q e (S),Q e (T)}. Moreover, 
it coincides with the quantum capacity on ideal channels: 
Qe(idn) = <3(id n ) = log 2 n, and it vanishes whenever 
TO is completely positive. In particular, Qq(T) = 0, 
whenever T is separable in the sense that it can be de- 
composed as T = PM into a measurement M and a 
subsequent preparation P based on the measurement re- 
sults. This follows immediately from the observation that 
on classical systems transposition is the identity. Then 
P© = QPQ is a channel, and so is MPQ. We note that 
Qe is also closely related to the entanglement quantity 
log 2 || (id (8) 6)(yo)||i, i.e., the logarithm of the trace norm 
of the partial transpose of the density operator, which 
enjoys analogous properties. 



III. QUANTUM GAUSSIAN STATES 
A. Canonical Variables and Gaussian states 

In this Section we recapitulate some results from [[l0| , 
||, @ for the convenience of the reader. Our approach 
to quantum Gaussian states is based on the characteris- 
tic function of the state which closely parallels classical 
probability ]i~l] ], p2| , and is perhaps the simplest and 
most transparent analytically. An alternative approach 
can be based on the Wigner "distribution function" p5] . 

Let QjiPj be the canonical observables satisfying the 
Heisenberg CCR 

[Qj,Pk] = iSjkhl, [qj,q k ] = 0, [p 3 ,Pk] = 0. 
We introduce the column vector of operators 

R = [qi,Pi, ■ 



lT 



the real column 2s-vector z — [xi,yi, . . . ,x s ,y s } T , and 
the unitary operators in TL 



V(z) = exp i^2(xjqj + yjpj) 
i=i 

= exp i R T z. 



(3.11) 



These "Weyl-operators" satisfy the Weyl-Segal CCR 

V(z)V(z') = exp[^A(z, z')]V(z + z% (3.12) 

where 



A(z, z') = nJ2( x 'jyj - xjy'j) 



(3.13) 



is the canonical symplectic form. The space Z of real 
2s-vectors equipped with the form A(z, z') is what one 
calls a symplectic vector space. We denote by 



% 
-h 



h 

-h o 



(3.14) 



the (2s) x (2s)-skew-symmetric commutation matrix of 
components of the vector R, so that 

A(z,z') = -z T Az'. 

Most of the results below are valid for the case where 
the commutation matrix is an arbitrary (nondegenerate) 
skew- symm etric matrix, not necessarily of the canonical 
form ( |3.14| ). 

A density operator p has finite second moments if 
Tr(pgJ) < oo and Tr(pp^) < oo for all j. In this case 
one can define the vector mean and the correlation ma- 
trix a by the formulas 



m = TrpR ; a — — A = Tr(i? — m)p(R 



(3.15) 



The mean can be an arbitrary real vector. The correla- 
tion matrix a is real and symmetric. A given a is the 
correlation matrix of some state if and only if it satisfies 
the matrix uncertainty relation 



-A > 0. 

2 ~ 



(3.16) 



We denote by E (m, a) the set of states with fixed mean 
m and the correlation function a. The density operator p 
is called Gaussian, if its quantum characteristic function 
(f>(z) = TrpV(z) has the form 



(z) = exp \ i m z 



1 



T 

z az 



(3.17) 



where m is a column (2s)-vector and a is a real symmetric 
(2s) x (2s)-matrix. One then can show that m is indeed 
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the mean, and a is the correlation matrix, and (3.17) 
defines the unique Gaussian state in E (m, a) . In what 
follows we will be interested mainly in the case to = 0. 

The correlation matrix a describes a quadratic form 
rather than an operator. Therefore its eigenvalues have 
no intrinsic significance, and depend on the choice of ba- 
sis in Z. On the other hand, the operator a defined 
by z T az — A(z, az) has a basis free meaning. In ma- 
trix notation it is a — A _1 a. This operator is always 
diagonalizable, and its eigenvalues come in pairs ±ijj- 
Diagonalizing this operator is essentially the same as the 
normal mode decomposition of the phase space, when the 
form z T az is considered as the Hamiltonian function of 
a system of oscillators. It leads to a decomposition of the 
phase space into two-dimensional subspaces, such that 
on the j th subspace we have (in some new canonical vari- 
ables qj,pj) 



7j 
7j 



1 

-1 



(3.18) 



and all terms between different blocks vanish. The ma- 
trix uncertainty relation now requires jj > 1/2, in which 
equality holds iff pj is the pure (minimum-uncertainty) 
state. Hence a general Gaussian state p is pure if and 
only if all jj — 1/2, or 



(A-M 2 = - \l, 
in which case E (to, a) reduces to a single point. 



(3.19) 



B. Gauge-invariant states 

We shall be interested in the particular subclass of 
Gaussian states most familiar in quantum optics, namely, 
the states having a P-representation 



9= I \0({\m(d 28 



(3.20) 



where ^jv(d 2s £) is the complex Gaussian probability 
measure with zero mean and the correlation matrix N. 
(see e.g. §§, Sec. V, 5. II). Here £ G C s , |<) are the co- 
herent vectors in 7i, a\Q — CIC)> N is positive Hermitian 
matrix such that 



N = Tr {a pa)) 



(3.21) 



(we use here vector notations, where a — [a\, . . . , a s ] T is 
a column vector and a) — [a\ , . . . , a\] is a row vector) 

and a J = 711 +i Po)- 

These states respect the natural complex structure in 
the sense that they are invariant under the gauge trans- 
formations a — > a exp(iip) . As shown in |l3| ] , the quantum 
correlation matrix of such states is 



a = h 



ReN + I/2 -ImJV 
ImN RcN + 1/2 



With Pauli matrices I%,o y , the real 2s x 2s— matrices of 
such form can be rewritten as complex s x s— matrices, 
by using the correspondence 



A -B 
B A 



= I 2 A + o-yB <-> A + iB, 
which is an algebraic isomorphism. Obviously, 

= Sp{A + iB), 



A -B 
B A 



where by "Sp" we denote the trace of matrices, as op- 
posed to the trace of Hilbert space operators, which is 
denoted by "Tr" . By using this correspondence, we have 



%(N + I/2), 



A 



-m, 



and 



i(N + I/2). 



(3.22) 



(3.23) 



For the case of one degree of freedom we shall be inter- 
ested in the last Section, N is just a nonnegative number 
and p is an elementary Gaussian state with the charac- 
teristic function 



(z) = exp 



N ■ 



(3.24) 



where we put \z\ 2 = (x 2 +y 2 )- This state has correlation 
matrix of the form ( |3.1g| ) in the initial variables q,p, with 
7 = N + 1/2, and is just the temperature state of the 
harmonic oscillator 



9i 



1 



7 + l/2^V 7 +l/2 



7-1/2V 



\n) (n\ 



(3.25) 



in the number basis \n), with the mean photon number 
N. 



C. Computation of entropy 

To compute the von Neumann entropy of a general 
Gaussian state one can use the normal mode decomposi- 
tion. For a single mode, the density operator pj with the 
correlation matrix (3.18), setting 7., = 7 for convenience, 



is unitarily equivalent to the state ( 3.25 ). From this one 
readily gets the von Neumann entropy H(py) by summa- 
tion of the geometric series, and for general Gaussian p 
by summing over normal modes. 

To write the result in compact form, one introduces 
the function 

g(x) = (x + l)log(x + l)-xlogx, x>0 (3.26) 
.9(0) = 0. 

Then 
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7J| 



(3.27) 



where jj runs over all eigenvalue pairs Hjj of A a. 

One can also write this more compactly, using the fol- 
lowing notations, which we will also use in the sequel. 
For any diagonalizable matrix M = Sdia,g(jrij)S~ 1 , we 
put abs(Af) = S'diag(|mj|)S'~ 1 , analogously for other 
conti nuou s functions on the complex plane. Then equa- 
tion ( 3.27 ) can be written as Jl3| 



H(P) = ^Sp g (abstA^a) - *- 



(3.28) 



For gauge-invariant state, by using ( 3.23 ), this reduces 
to the well-known formula 

H(p) = Sp g{N). 



D. Schmidt Decomposition and Purification 

Forming a composite systems out of two systems de- 
scribed by CCR-relations is very simple: one just joins 
the two sets of canonical operators, making operators 
belonging to different systems commute. The symplec- 
tic space of the composite system is a direct sum Zx2 = 
Z\®Z2, which means that elements of this space are pairs 
(zi, z 2 ) with components Zi G Zi. In terms of Weyl oper- 
ators one can write V12 (21,22) = Vi(2i)®V^(.Z2). By def- 
inition, the symplectic matrix A12 is block diagonal with 
respect to the decomposition Z = Z\ ® Z2. However, the 
correlation matrix a±2 is block diagonal if and only if the 
state is a product. The restriction of a bipartite Gaussian 
state p to the first factor is determined by the expecta- 
tions of the Weyl o perators Vj (21) 8 1 = 1^2(21, 0), hence 
according to (3.17), by the the correlation matrix ax with 



zfaizi = (21, 0) 0:12(21, 0), which is just the first diag- 
onal block in the block matrix decomposition 



01X2 



a x (3 
3 T 



o: 2 



A12 = 



Ai 
A 2 



(3.29) 



As in the case of bipartite systems with finite dimen- 
sional Hilbert spaces there is a canonical form for pure 
states of the composite system, the Schmidt decomposi- 
tion. Like the diagonalization of a one-site density oper- 
ator, it can be carried out for Gaussian states at the level 
of correlation matrices. By writing out equation (3. IE) 
in block matrix form, we find in particular that 



(Ar 1 ai)(A r 1 /3) = (Ar 1 /3)(A 2 - 1 a 2 ). 



(3.30) 



Thus (Aj" 1 /^) maps eigenvectors of (A 2 ~ 1 «2) into eigen- 
vectors of (A^j~ 1 ai), with the same eigenvalue. Hence the 
spectra of the restrictions are synchronized much in the 
same way as in the finite dimensional case, and all the 



matrices ax , a.2 , 3 can be diagonalized simultaneously by 
a suitable choice of canonical coordinates. Evaluating 



also the diagonal part of Equation (3.19), one gets an 
equation for 0, so that finally ai2 is decomposed into 
blocks corresponding to (a) pure components belonging 
to only one subsystem, and not correlated with the other, 
and (b) b locks of a standard form, which can be w ritte n 
like ( j3~29l ) with ax = a 2 = a, A x = A 2 = A from (pU8| ), 
and 



r 



1 

-1 



(3.31) 



The purification of a general Gaussian state can easily 
be read off from this, by constructing such a standard 
form for every normal mode. In order to write in 
operator form without explicit reference to the normal 
mode decomposition, it is most convenient to perform 
an appropriate reflection in the space Z 2 , by which 
becomes purely off-diagonal. Then we can choose |27j 
Ai = A = — A2 and a 2 = ax = a, resulting in 



3 = -B T = AV~(A- 1 a) 2 -//4. 



(3.32) 



This also covers cases with 3 = for some modes, where, 
strictly speaking no purification would have been neces- 
sary. We thus have 



A^a^ 



A^a 



V-(A-ia)2-//4 



vMA-ia) 2 -1/4 -A- X a 



(3.33) 



In the gauge-invariant case, we can use the correspon- 
dence 



Ajax2 



i(N + I/2) VWTN 



VW+N -i(N + I/2) 
following from ( [3.22 ) . 



(3.34) 



IV. LINEAR BOSONIC CHANNELS 

A. Basic Properties 

The characteristic property of the channels considered 
in this paper is their simple description in terms of phase 
space structures. The key feature is that Weyl operators 
go into Weyl operators, up to a factor. That is, the 
channel map in the Heisenberg picture is of the form 



T*(V'(z')) = V(K T z')f(z'), 



(4.35) 



where K : Z — > Z' is a linear map between phase spaces 
with symplectic forms A and A', respectively, and f(z') 
is a scalar factor satisfying certain positive definiteness 
condition to be discussed later. Because of the linearity 
of K, such channels are called linear Bosonic channels 



G 



|L5||, and if, in addition, the factor / is Gaussian, T will 
be called a Gaussian chan nel. In terms of characteristic 
functions, Equation flf.35| ) can be written as 



0V) - <P(K T z')f{z') 



(4.36) 



where <f> and <j)' are the characteristic functions of input 
state p and output state T[p], respectively. 
We will make use of following key properties: 

(A) The dual of a linear Bosonic channel transforms 
any polynomial in the operators R' into a polynomial in 
the R of the same order, provided the function / has 
derivatives of sufficiently high order. This property fol- 
lows from th e de finition of moments by differentiating 
the relation ( 4.35f ) at the point z' — 0. 

(B) A Gaussian channel transforms Gaussian states 
into Gaussian states. This follows from the definition of 



Gaussian state and the relation (4.36). 

(C) Linear Bosonic channels are covariant with re- 
spect to phase space translations. That is if p z — 
V{~A- 1 z)pV(-A- 1 z)* is a shift of p by z, T[p] is simi- 
larly shifted by Kz. 

There is a dramatic difference in the capacities of a 
Gaussian channel for classical as opposed to quantum in- 
formation. Classical information can be coded by using 
phase space translates of a fixed state as signal states, so 
the output signals will also be phase space translates of 
each other. Then no matter how much noise the channel 
may add, if we take the spacing of the input signals suf- 
ficiently large, the output states will also be sufficiently 
widely spaced to be distinguishable with near certainty. 
Therefore the unconstrained classical capacity is infinite. 
The same would be true, of course, for a purely classical 
channel with Gaussian noise. The classical capacity of 
such channels becomes an interesting quantity, however, 
when the "input power" is taken to be constrained by a 
fixed value, which we must take as one of the parame- 
ters defining the channel. Then arbitrarily wide spacing 
of input signals is no longer an alternative, because an 
intrinsic scale for this spacing has been introduced. 

The remarkable fact of quantum information on Gaus- 
sian channels is that such an intrinsic scale is already 
there: it is given by h. As we will show, the quantum in- 
formation capacity is typically bounded even without an 
energy constraint. Loosely speaking, although we send 
arbitrarily many well distinguishable quantum signals 
through the channel, coherence in the form of commu- 
tator relations is usually lost. Surprisingly, in spite of 
the infinite classical capacity, the capacity for quantum 
information may be zero, which means that even joining 
arbitrarily many parallel channels with poor coherence 
properties is not good enough for sending a single qubit. 
This phenomenon will be explained in some detail in Sec- 
tion V. 

The choice of the scalar function f(z') is crucial for the 
quantum transmission properties of the channel. Nor- 
malization of T requires that /(0) = 1, and it is clear 



Beyond that, it is not so easy to see which choices of / 
are compatible with the complete positivity. If / decays 
rapidly, T* maps most operators to operators near the 
identity, which means that there is very much noise. On 
the other hand, there will be a lower limit to the noise, 
depending on the linear transformation K . Only when 
if is a symplectic linear map and T is reversible, the 
choice f(z) = 1 is possible. Otherwise, there is some 
unavoidable noise. 

There are two basic approaches to the determination of 
the admissible functions /. The first is the familiar con- 
structive approach already used in Section II, based on 
coupling the system to an environment, a unitary evo- 
lution and subsequent reduction to a subsystem, with 
all of these operations in their linear Bosonic/Gaussian 
form. Basically this reduces the problem to linear trans- 
formations of systems of canonical operators. This will 
be described in Subsection B, and used for the calcula- 
tion of entropy exchange in Subsection C. Alternatively, 
one can describe the admissible functions / by a twisted 
positive definiteness condition, and this will be used for 
evaluating the bound Cq(T) in Subsection D. 



B. Bosonic channels via transforming canonical 
operators 

Let R,Re be vectors of canonical observables in 
H.,H.e, with the commutation matrices A, A^;. Consider 
the linear transformation 



R' = KR + KeRe 



(4.37) 



where K, Ke are real matrices (to simplicfy notations 
we write R, Re instead of iJ® /e, 7(8 Re etc.) Then the 
commutation matrix and t he correlation with respect to 
R' are computed via ( pM5| ) with m = 0, namely 



a 



I J U'T 



-A' = TrR'p'R 
2 p 



We apply this to the special case p' — p <g> pe, where 
Pe and p are density operators in H.e and H. with the 
correl ation matrices cxe and a, respectively. Then using 
(4.37), we obtain 



A' = KAK T + K E A E K^ 
a' = KaK T + K E a E K^. 



(4.38) 



that |/(z')l < 1 for all z' , from taking norms in (4.35) 



Of course, the operators R' need not form a complete set 
of observables in H <8> He, but in any case a' is the cor- 
relation matrix of a system containing just the canonical 
variables R' , and it is this state which we will consider 
as the output state of the channel. 

For fixed state pe (state of the "environment") the 
channel transformation taking the input state p to the 
output p 1 is described most easily in terms of character- 
istic functions: 
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4>'{z') = cp{K T z')^E (K T E z') . (4.39) 
We can write this as a linear Bosonic channel in the form 



fl4.36| ) with 

/CO 



Trp E V E (K E z') 



(4.40) 



Thus the factor / is expressed in terms of the charac- 
teristic function of the initial state of the environment. 
Obviously, the channel is Gaussian if and only if this state 
is Gaussian. 

If we want to get the state of the environment after 
the channel interaction, as required in the definition of 
exch ange e ntropy, we have to supplement the linear equa- 
tion (4.37) by a similar equation specifying the environ- 



ment variables R' E after the interaction: 

R' = KR + K E R E , 
R'e = LR + LeRe. 

Assuming that Z = Z' and A' = A, one can always 
choose L,Le such that the combined transformation is 
canonical, i.e., preserves the commutation matrix 



A 
A E 



Then the channel Te 
relation 



P 



p' E can be defined by the 



T* E [Ve{ze)] = V{L T z E ) ■ 4>e {L T E z E ) , 
and is thus also linear Bosonic. 

C. Maximization of mutual information 

The estimate for the entanglement assisted classical 
capacity suggested by |ll| is t he maximum of the quan- 
tum mutual information ( |2 .3| ) over all states satisfying 
an appropriate energy constraint. Evaluating this max- 
imum becomes possible by the following result :Q Let T 
be a Gaussian channel. The maximum of the mutual in- 
formation I(p) over the set of states £ (to, a) with given 
first and second moments is achieved on the Gaussian 
state. 

Proof (sketch). By purification (if necessary), we can 
always assume that pe is pure Gaussian. Then we can 
write 

I(p) = H(p) + H(T[p])-H(T E [p}). 

Let po be the unique Gaussian state in £ (to, a) . For 
simplicity we assume here that po is nondegenerate. The 



1 The proof of this theorem was stimulated by a question 
posed to one of the authors (A.H.) by P. W. Shor. 



general case can be reduced to this by separating the pure 
component in the tensor product decomposition of po- 
The function I(p) is concave and its directional derivative 
at the point po is (cf. [fl8|| ) 

VxHpo) = TrX(lnp + I) + Tr T[X}(hxT\p ] + I) 
-TrT E [X](\nT E [p ] + I). 

By using dual maps this can be modified to 

V x /(po) = TrX{lnp + T*[lnT[p }] 

-T E []nT E [p ]]+l}. (4.41) 

Now by property (B) of Gaussian channels, the opera- 
tors pa,T[po],T E [po] are (nondegenerate) Gaussian den- 
sity operators, hence their logarithms are quadratic poly- 
nomials in the corresponding canonical variables (see Ap- 
pendix in |13j). B y property (A) the expression in curly 
brackets in (4.41) is again a quadratic polynomial in R, 
that is a linear combination of the constraint op erat ors 
in £ (m, a). Therefore, the sufficient condition (1.3) in 
the Appendix is fulfilled and I(p) achieves its maximum 
at the point po € £ (to, a). 

This theorem implies that the maximum of I(p) over a 
set of density operators defined by arbitrary constraints 
on the first and second moments is also achieved on a 
Gaussian density operator. In particular, for an arbi- 
trary quadratic Hamiltonian H the maximum of I(p) over 
states with constrained mean energy TrpH is achieved on 
a Gaussian state. The energy constraint is linear in terms 
of the correlation matrix: 



Sp(ea) < N, 

where e is the diagonal energy matrix (see |l3|]). 

When p and T are Gaussian, the quantities H(p), 
H(T[p]),H(p,T) and I(p,T),J(p T) can in principle 
be computed by using formulas (|3~28| ), (^38|) , ( p3| ). 
Namely, H(T[p\) is given b y for mula (3.28) with a re- 
placed by a 1 computed via ( 1.35 ) , and 



H(p,T) = is P3 (abs(Ar 2 



Q, 



2 ) 



where 



a 

I3 T K T 



a 



13 = Ay/-(A-h 



J/4 



is computed by inserting ( 1.37 ) into 



l 12 



A' 12 =Ty(R',R 2 )p(R',R 2 ) 



where i?2 are the (unchanged) canonical observables of 
the reference system. 

Alternatively, the entropy exchange can be calculated 
as the output entropy H(T E [p\) if an explicit description 
of Te is available. We shall demonstrate this method in 
the example of one-mode channels in the Appendix. 
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D. Norms of Gaussian Transformations 

The transposition operation on a Bosonic system can 
be realized as the time reversal operation, i.e., the opera- 
tion reversing the signs of all momentum operators, while 
leaving the position operators unchanged. Obviously, the 
dual T* then takes Weyl operators into Weyl operators. 
So transposition is just like a linear Bosonic chann el, al - 
beit without the scalar factor f(z') in Equation ( 4.36 ). 
It is this factor which makes the difference between pos- 
itivity and complete positivity, and also enters the norm 
||T|| c b. In this Subsection we will provide general criteria 
for deciding complete positivity and computing the norm 
of general linear Bosonic transformations. 

These are by definitio n the operators T acting on Weyl 
operators according to ( [4.35| ) where f(z') is a scalar fac- 
tor. We will assume for simplicity (and in view of the 
applications in the following sections) that the antisym- 
metric form 

A"( Zl ,z 2 ) = A'(2i, z 2 ) - A(K T Zll K T z 2 ) (4.42) 

is non-degenerate. This makes the space Z' with the 
form A" into a phase space in its own right. With 
the introduction of suitable canonical coordinates it be- 
comes isomorphic to (Z, A), so there exists an invert- 
ible linear operator A : Z — » Z such that A"(z\,z 2 ) = 
A(A- 1 z 1 ,A- x z 2 ). 

If / is continuous and has sufficient decay properties 
(which will be satisfied in our applications), there is a 
unique trace class operator p determined by the equation 



Tv(pV(z)) = f(Az). 



(4.43) 



Then T is completely positive if and only if p is a pos- 
itive trace class operator. This is a standard result in 
the theory of quasi-free maps on CCR-algebras |J. It is 
proved by showing that both properties are equivalent to 
a "twisted positive definiteness condition" , namely the 
positive definiteness of all matrices of the form 

M rs = f(z r ~z s ) exp(-^A'(z r ,z s ) + ^A(K T z r ,K T z s )), 

where z\ , . . . , z n are an arbitrary choice of n phase space 
points. 

If p is a non-positive hermitian trace class operator, 
it has a unique decomposition into positive and neg- 
ative part: p = p + — p_ such that p± > 0, and 
p+P- = 0. Then \p\ — p + + p_ and the trace norm 
is = Tr(p + ) + Tr(p-). Inserting p± into Equa- 



tion (4.43) instead of p, w e get two functions f± on phase 
space and from Equation (4.35) two linear Bosonic trans- 



formations T± with T = T + — T_. By the criterion just 
proved, T+ and T_ are completely positive. Hence 

||T|| c b < HT+IU + ||r_|| cb = ||T + (1)|| + ||T_(1)|| 

= /+(0) + /-(0) = Tr( P+ ) + Tr(p_) = ||p|| a (4.44) 

If the factor / is a Gaussian, i.e., 



f(z) = cxp(--z I3z) 



(4.45) 



for some positive definite matrix /3, we can go one step 
further. In this case we may decompose (3 into normal 
modes with respect to A", which decomposes T into a 
tensor product of one- mode Gaussian transformations Tg, 
for each of which ||X^|| c b may be computed separately by 
the above method. This amounts to com puting the trace 
norm of the operator p 7 given by ( 3.25 ) with arbitrary 
positive 7. The absolute value of p 1 is obtained by taking 
absolute values of all the eigenvalues, which still makes 
||p 7 ||i a geometric series: 



1 



IP7II1 = 



7+1/2 



E 



7-1/2 



7+1/2 



1 



= max{l, — }. (4.46) 
2 7 



This is all the information we need for the estimates of 
quantum capacity in the following Section. 



V. THE CASE OF ONE MODE 



A. Attenuation/amplification channel with classical 



The channel we consider in this Section combines at- 
tenuation/amplification jl4| with additive classical noise 
ifLsf . It can also be described as the most general one- 
mode gauge invariant channel, or in quantum optics ter- 
minology, the most general one-mode channel not involv- 
ing squeezing. Channels of this type were also used in 
fl28| as the basis for an analysis of the classical limit of 
quantum mechanics. 

Let us consider the CCR with one degree of freedom 
a = -^=^(<Z + ip), and let do be another mode in the 
Hilbert space Tto = He of an "environment". Let the 
environment be initially in the vacuum state, i.e., in the 
state with the characteristic function ( |3.24 ) with N = 
0. Let £ be a complex random variable with zero mean 
and variance N c describing additive classical noise in the 
channel. The linear attenuator with coefficient k < 1 and 
the noise iV c is described by the transformation 



a' = ka + y/l — k 2 ao + £ 

in the Heisenberg picture. Similarly, the linear amplifier 
with coefficient k > 1 is described by the transformation 



a = ka+ \Jk 2 - laj + £. 

It follows that the corresponding transformations T[p] of 
states in the Schrodinger picture both have the charac- 
teristic function 



TrT[p}V(z) = TrpV(kz) x 

h 



x exp 



(\k 2 -l\/2 + N c )\zf 



(5.47) 
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Let the input state p of the system b e the elementary 
Gaussian with characteristic function ( 3.24 ). Then the 
entropy of p is H(p) — g(N). From ( |5.47 ) we find that 
the output state T[p] is again elementary Gaussian with 
N replaced by 



where 



N' = k 2 N + Nq, 



N[) = max{0, {k 2 - 1)} + iV c 



is the value of the output mean photon number corre- 
sponding to the input vacuum state. Then 



H(T[p])=g(N>). 



(5.48) 



Now we calculate the exchange entropy H(p,T). The 
(pure) input state p\2 of the extende d syst em Hi ® Hi 
is characterized by the 2 x 2— matrix ( 3.34 ). The action 
of the extended channel (T <S> id) transforms this matrix 
into 



i(N>+l) 



_ ky/N(N+l) 

1) -i(iV+|) 



From formula ( ^27| ) we deduce H(p,T) = g(\Xi\ - 5) + 
5(1^21 — \): where Ai, A2 are the eigenvalues of the com- 
plex matrix in the right-hand side. Solving the charac- 
teristic equation we obtain 



{(N' -N)±D), 



(5.49) 



where D = J (N + N' + 1) - 4k 2 N(N + 1). Hence 



H(p,T) 



D + N' - N - 1 



9 



(5.50) 
D - N' + N - V 



Now using the theorem of Section 5, we can calculate 
the quantity 

C e (T) = I(p, T) = H{p) + H(T[p}) - H(p, T) 

as a function of the parameters N,k,N c , and try to com- 
pare it with the one-shot unassisted class ical capacity of 
the channel Ci(T) given by expression (2.4) where the 
maximum is taken over all probability distributions {pt} 
and the collections of density operators {pt}, satisfying 
the power constraint ^^^piTrpia^a < N. It is quite plau- 
sible, but not yet proven that this maximum is achieved 
on coherent states with the Gaussian probability density 
p{z) = (irN)^ 1 exp (— \z\ 2 /N), giving the value 



C 1 (T) = g(N')-g(N^ 



The ratio 



G = 



(5.51) 



then gives at least an upper bound for the gain of using 
entanglement-assisted versus unassisted classical capac- 
ity. In particular, when the signal mean photon number 

>0, 



N tends to zero while Nq 



C X {T) ~ Nk 2 log 



K 



C e {T)~ -N log N/W + l), 

and G tends to infinity as — log N. 

The plots of G as function of k for N c = 0, and 
as a function of N c for k — 1 are given in Figure 1 
and Figure 2, respecitively. The behavior of the en- 
tropies H(T[p\), H(p,T) as functions of k for N c = is 
clear from Figure 3. For all N the coherent information 
H(T[p\) - H(p, T) turns out to be positive for k > l/y/2 
and negative otherwise. It tends to —H(p) for k — > 0, 
is equal to H(p) for k = 1, and quickly tends to zero as 
k — > 00 (see Figure 4). 



B. Estimating the quantum capacity 

Going back to the upper bound for quantum ca pacity 
in Section IV, we see that T is given by equation ( 4.35 ) 
with Kz = kz and 



f(z) = exp(- 



(|fc 2 -l|/2 + iV c ) 



M 2 ). 



Then A" = (1 — fc 2 )A, and the operator A mapping the 
symplectic form A" to the standard form A is multipli- 
cation by \J\k 2 — 1|, combined for k > 1 with a mirror 
reflection to change the sign. This leaves 



f{Az) 



exp 



(|fc 2 -i|/2 + Ag 

2|fe 2 -l| 



\z\' 



(5.52) 



i.e., p = p 1 with equations (3.25) and ( 4.43| ), where 
7 = 1/2 + N c /\k 2 - 1|. This is the verification of the 
complete positivity of T by the methods of the above 
section. Of course, this is strictly speaking unnecessary, 
because T was constructed explicitly as a completely pos- 
itive operator in terms of its dilation in Subection IV. A. 

But let us now consider TO. It is also a Bosonic linear 
transformation, in which O only has the effect of chang- 
ing the sign of the symplectic form, without changing /. 
Thus A" = (1 + k 2 )A, and 



f(Az) = exp 



(|fc 2 -l|/2 + jV c ) 
2\k 2 + l\ 



4 



which see ms like a rather minor change over Equa- 
tion (5.52). However, we now get p = p 1 with 7 = 
(|fc 2 -l|/2 + 7V c )/(fc 2 + l) which is not necessarily > 1/2, 
so TO is not necessaril y com pletely positive. Taking the 
logarithm of Equation ( 4.46 ) we get 



Q e (T) <max{0, 
log 2 (fc 2 



l)-log 2 (|fc 2 -l| + 27V c )}. (5.53) 
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In particular, for 7 > 1/2, i.e., for N c > (\k 2 + 1| - |fc 2 - 
l|)/2 = max{l,fc 2 }, the capacities Q®(T), and hence 
g e (T) and Q(T) all vanish. 

This upper bound on quantum capacity is interest- 
ing to compare with the quantity Qq(T) = sup J(p 7 T), 
where J{p 1 T) = H (T[p\) — H (p,T) , and the supremum is 
taken over all Gaussian input states. Since the coherent 
information 



J(p,T) = g(N') - g ( 



D + N'-N-l 



D-N' + N-l 



(5.54) 



increases with the input power N, we obtain 

Q G {T) = lim J{p,T) (5.55) 



N- 



= logfc 2 -log|fc 2 -l|-. 9 (A c /|fc 2 -l|), 

which is in a good agreement with the upper bound 
( 5.53 ) (see Figure 4). 



If (1 — t)po + ta > for small negative t, then we say 
that the direction apo is inner. In that case ([O]) takes 
the form 



VxF(po) = 0. 



(1.2) 



If po is nondegenerate, then the direction apo is inner for 
arbitrary pure a in the range of yfpo, and the neces sary 
and sufficient condition for the minimum is that (1.2) 
holds for arbitrary such a. 

Let Ai, i = 1, . . . , r be a collection of selfadjoint con- 
straint operators. Assume that for some real constants 
A* 



VxF{p ) =TrX J^XiAi. 



(1.3) 
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APPENDIX A: MINIMIZING CONVEX 
FUNCTION OF A DENSITY OPERATOR. 

There is a useful lemma in classical information the- 
ory which gives necessary and sufficient conditions for 
the global minimum of a convex function of probability 
distributions in terms of the first partial derivatives. The 
lemma is based on general Kuhn- Tucker conditions and 
can be generalized to functions depending on density op- 
erators rather than probability distributions. 

Let F be a convex function on the set of density op- 
erators S, and pa a density operator. In order F to 
achieve minimum on po, it is necessary and sufficient 
that for arbitrary density operator a the convex function 
F((l — t)po + ta) of the real variable t achieves minimum 
at t — 0. For this, it is necessary and sufficient that 



Let us consider CCR with one degree of freedom de- 



VxF(po) = 



dt 



F((l-t) Po +ta)>0, (1.1) 



where X = a — po, and V xF(po) is the directional deriva- 
tive of F in the direction A, assuming that the deriva- 
tives exist. If a = J^iPi °^ tnen ^xF(p ) = J2iPi 
V XiF(po), where X, = at— po. Therefore it is necessary 
and sufficient that ( O ) holds for pure a. 



2/1 



scribed by one mode annihilation operator a 
ip), and consider the transformation 

a = a + £, 



where £ is a complex random variable with zero mean 
and variance N r . This is a transformation of the type 



(4.37) with Ag — 0, which describes quantum mode in 
classical Gaussian environment. The action of the dual 
channel is 

T*[f(a,J)}= [ f{a + z,(a + z)^)^ Nc (d 2 z), 



where z = -j=(x + iy) is now complex variable, and 

/Ltjv c (ci 2 z) is complex Gaussian probability measure with 
zero mean and variance N c , while the channel itself can 
be described by the formula 



T\p] 



D(z)pD(z)* mc (d 2 z), 



(2.1) 



where D(z) = expi (za^ — za) is the displacement oper- 
ator. 

The entanglement-assisted classical capacity of the 
channel (|]l]) was first studied in |Q by using rather 
special way of purification and the computation of the 
entropy exchange. A general approach following the 
method of was described in Sections IV-V; here we 
give an alternative solution based on the computation of 
the environment entropy. 

For this we need to extend the environment to a quan- 
tum system in a pure state. Consider the environment 
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Hilbert space He — L 2 ([ijyJ with the vector |^o) given 
by the function identically equal to 1. The tensor prod- 
uct H ® He can be realized as the space L^(fj,j^ c ) of 
/iAr c -square integrable functions ip(z) with values in H. 
Define the unitary operator U in H (g) He by 



Thus 



Then 



while 



(UiP)(z) = D(z)i>(z). 
T{p} = TY HE U(p®\*o)(*o\)U*, 
TE[p]=Tr H U(p®\V )(y \)U*. 



This means that Te[p] is an integral operator in L 2 (/i/v c ) 
with the kernel 

K{z,z') = Tr D(z)p a D(z')* 

= cxp(i%z'z -(E+ l/2)\z - z'\ 2 ). 

Let us define unitary operators V(z%, z%) in L 2 (/^jv ( .) by 

V(z\,zi)il>{z)= tf)(z + z 2 ) x 

x exp 



im(z + y) - J^m(z + y) 



The operators V(zi, z<x) satisfy Weyl-Segal CCR with two 
degrees of freedom with respect to the symplectic form 

A((«i, z 2 ), (z[,z' 2 )) = R {z[z 2 ~ z x z 2 ) . 

Passing over to the real variables x, y one finds the cor- 
responding commutation matrix 



0-10 
-1 
10 
10 



The characteristic function of the operator Te[po] is 
TyTe[po}V(z 1 ,z 2 ) = / V(z 1 ,z 2 )K(z,z)\- =z p Nc (d 2 z), 



where V(zi, z 2 ) acts on K as a function of the argument 
z. Evaluating the Gaussian integral, we obtain that it is 
equal to 



exp 



D 2 

2N c % Zl z 2 + — 1 2 2 



(where now D = ^J{N C + l) 2 + 4N C N), which is Gaus- 
sian characteristic function with the correlation matrix 



N c 

iV c 

-N c 

N c 



N c 
-N c 
^ 

o eL 



A E a E= 2 





N c 
-N c 




-N S± 







-Nr. N, 



SL 



-JV C 





By using Pauli matrix a y , we can write it as 



-ia y N c 



I 

a y 



i&yN c <7 y N ^ 

-a y N c -ia y N c 



I 

a v 



hence the absolute values of the eigenvalues of A E a' E 
are the same as that of the matrix 



iN r 



N r iN r 



which coincide with (|5.49| ) in the case k = 1 
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FIG . 1. Gain of entanglement assistance. 
Gain (5.51) as a function of k with N c = 0. Parame- 
ter=input noise N. 
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Nc 



FIG . 2. Gain of entanglement assistance. 
Gain ( 5.51 ) as a function of N c with k = 1. Parame- 
ter=input noise A'". 




FIG. 3. E ntropies. 
outpu t entropy from ( |5.4g| ), exchange entropy from 
(E50|) with N c = 0. 
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FIG. 4. Bounds for Quan tum Capacity, N c 
J=coherent information ( 5.54 ) with N = .7; 



0. 



QG=Qg= bound maximized over Gaussians (5.55); 

QT=bound Qe from transposition ( 5.53 ); 

Z= zero at k = l/v2, common to all curves of type J. 




FIG. 5. Gaussian maximized coherent information Qc{T) 
as function of k and N c . The shaded area is the area, where 
Qe > 0. 



